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ABSTRACT 

We present results of Monte Carlo calculations for the Lya line transfer in an 
expanding dusty supershell, where Lya source is a well-localized star cluster in 
a starburst galaxy. The escape of Lya photons from such system is achieved by 
a number of back-scattering, and so a series of emission peaks are formed red- 
ward of the systemic redshift by back-scattering. However, majority of observed 
Lya emission from starbursts show singly-peaked asymmetric profiles. We find 
in this paper that, in order to form a singly-peaked Lya emission, dust should 
be distributed in the ionized bubble, as well as within the supershell of neutral 
hydrogen. We also find that the overall escape fraction of Lya photons is de- 
termined by the HI column density of the supershell, the expansion velocity of 
the supershell, and the spatial distribution of dust. However, the kinematic in- 
formation of the expanding supershell is preserved in the profile of Lya emission 
even when the supershell is dusty. Our results are potentially useful to fit the 
P-Cygni type Lya line profiles from starburst galaxies, either nearby galaxies or 
high-z Lyman break galaxies (LBGs). 

Subject headings: line: profiles - radiative transfer - method: numerical - galax- 
ies: starburst 



1. Introduction 

Lya is the most prominent line feature in the rest-frame ultraviolet (UV) spectra of 
starburst galaxies, and so Lya is often used as a redshift indicator. Recently a large number 
of star-forming galaxies have been spectroscopically observed, either by using the Lyman 
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break method (Steidel et al. 1996, 1999; Shapley et al. 2003) or by using gravitational lens 
(Franx et al. 1997; Frye et al. 2002). The detailed review was given by Taniguchi et al. (2003). 
The rest-frame UV spectra of those galaxies often show unique emission. In these cases, we 
usually identify this emission as Lya, and use them as a redshift indicator. Moreover, we 
propose that this unique feature has more sources of astrophysical information of those 
star-forming galaxies. 

Lya profiles of starburst galaxies can be classified into three types: P-Cygni type emis- 
sion, broad absorption, and symmetric emission. The P-Cygni type emission consists of 
absorption in red part and the asymmetric emission. However, the rest-frame UV continua 
of high- 2; galaxies are often too weak to be detected, and so we can only see the asymmetric 
emission lines. The fraction of symmetric emission is small when compared with the other 
two types, and thus we will not go into detail of the symmetric type. 

Kunth et al. (1998) and Shapley et al. (2003) showed that the kinematics of ambient 
material is crucial to determine whether emergent Lya has either broad absorption or P- 
Cygni type emission. P-Cygni type Lya emission can be seen only in those galaxies whose 
interstellar absorption lines are blueshifted with respect to either Lya emission or stellar 
atmospheric lines. It seems to be a general consensus that outflowing motion is caused 
by either galactic superwind or galactic supershell (Lee & Ahn 1998; Heckman et al. 1998; 
Taniguchi & Shioya 2001; Ahn & Lee 2002; Ahn, Lee, & Lee 2003a,b; Silich et al. 2003; 
Spinrad 2003). The formation of asymmetric Lya by outflows was discussed qualitatively 
by Tenorio-Tagle et al. (1999). 

We study in this paper the formation of P-Cygni type Lya. Lya radiative transfer is 
mainly achieved by back-scattering processes in outflowing media. We have studied the role 
of back-scattering processes in the formation of Lya line profile (Ahn, Lee, & Lee 2003a,b). 
For the case of dust-free supershell, we found that a series of peaks appear in the red part 
of Lya. In that paper we emphasized that the kinematic motion is imprinted on both the 
width of each peaks and the velocity difference between peaks. However, majority of the 
observed Lya emission lines are dominantly singly peaked. Possibly we can attribute this 
discrepancy to the previously neglected dust extinction in the radiative transfer. We study 
in this paper the Lya line formation mechanism by interplay among the Lya resonance 
scattering in optically thick supershell of neutral hydrogen, the spatial extension of dust in 
the HII bubble, and the kinematics of the supershell with respect to Lya sources at the 
center of the supershell. 
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2. Model 

Galactic supcrshclls are very well-known structures in the nearby starburst galaxies 
(Marlowe et al. 1995; Martin 1998; Korthes & Kerton 2002). Their origin is known to be 
multiple explosions of the supernovae in active star-forming regions. In this paper we adopt a 
model galaxy of high- 2; starburst galaxies in which Lya sources are presumed to be located at 
a restricted part of the galaxy and surrounded by a dusty galactic supershell. The supershell 
is assumed to be made of uniform medium of neutral hydrogen and expanding in a bulk 
manner. We assume that the interior of the supershell is fully ionized, which is vacuum 
in the sense of Lya scattering. The model in this study is very simple, and it has several 
model parameters. However, the main purpose of this paper is to emphasize the role of 
spatial distribution of dust in destroying the secondary and higher peaks. We will itemize 
and clarify the assumptions in this study as follows. 

1. The starbursting region is assumed to be localized in a region of the model galaxy. 
This assumption is based on the fact that almost all the observed Lya emission lines 
are asymmetric. If the galaxy were rotating and its star-forming regions are scattered 
along spiral arms, we might see symmetric Lya emission lines. Hence we can consider 
that the star-forming regions are locahzed in a rather narrow region in the galaxies. 
The HST images of about half of the LBGs show the amorphous shape of those high-^ 
galaxies, while the other half show the nuclear starburst (Calzetti & Giavalisco 2001). 
Here we assume that even in the amorphous galaxies, only the brightest star-forming 
region is a dominant source of Lya photons. 

2. The ionized inner part of the supershell can also be a source of Lya photons. However, 
the fraction of photons, generated in the ionized part of the supershell, is smaller than 
that of the stellar Lya photons (Tenorio-Tagle et al. 1999). Moreover, there is no 
kinematic difference between the photons back-scattered by the neutral supershell and 
the photons emitted by ionized supershell. Hence we will neglect photons emitted from 
the ionized inner part of the supershell. 

3. We adopt a typical width of input Lya profile cri^ya = 50 km s~^ in this paper. Legrand 
(1997) observed B.a emission of a nearby HII galaxy Haro 2. The Ha line width is 
c"Ha = 50 km s^^ which is equivalent to 0"Lya = 50 km s~^. The input line width of 
Lya affects the width of emergent emission peaks. However, it ranges between 50 and 
115 km s~^ (Pettini et al. 2001), which justifies our choice. 

4. The continuum is neglected in this paper. According to Shapley et al. (2003), the Lya 
profiles for 38% of their sample of LBGs consist of a combination of both emission 
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and absorption. The equivalent width (EW) of broad absorption is 10 — 20A, which 
will be dumped in the red wing of the primary peak. However, those photons must be 
partially destroyed by dust extinction. Moreover, dust-free galaxies would have Lya 
equivalent widths of 50 — 200A (Chariot & Fall 1993), which is much larger than the 
absorption EW of the continuum. Hence we assume that the continuum will give a 
minor contribution. 

5. We also assume that the galactic supcrshcll fully covers the starburst region, and that 
Rmin = 0.9Rmax- Here Rmax IS the outer radius of the supershell and Rmin is the inner 
radius of the supershell. 

6. We also assume that the HI column density outside the supershell is negligible, because 
the supershell is assumed to be of galactic scale. We adopt the column density of neutral 
hydrogen A^hi = 10^^ — 10^°cm^^ in accordance with Voigt fittings done by Kunth et 
al. (1998). This column density corresponds to the line-center optical depth by Tq = 
2.27 (6/80 km s'^)" cm ]. Here b is the Doppler parameter corresponding 
to turbulence, which is typically 80 km s~^ for nearby HIl galaxies (Kunth et al. 1998). 
This gives us the Voigt parameter a = 7.60 x 10~^(80 km s~^/b). 

7. We also assumed for simplicity that the neutral hydrogen in the supershell is uniformly 
distributed. In fact, the supershell can either be patchy or have density gradient, which 
is also important for escape of Lyman limit photons. The geometry of the supershell is 
assumed to be spherical for simplicity. The supershell evolves more rapidly toward the 
steepest density gradient (Sihch & Tenorio-Tagle 1998), but it evolves to be of galactic 
scale. Thus we can assume that the supershell is spherical. However, more work for 
complicated geometry can be expected in future work. 

8. We adopt the typical expansion velocity of the supershell Vexp = 200 km s~^. Accord- 
ing to Kunth et al. (1998), the neutral gas in nearby starburst galaxies along the line 
of sight is being pushed by an expanding envelope around the HH region, outflowing at 
velocities close to 200 km s~^. Moreover, in the composite spectrum of high- 2; galax- 
ies, Shapley et al. (2003) measured an average blueshift for the strong low-ionization 
interstellar features of A^; = —150 ± 60 km s"-*^ with respect to the stellar systemic 
redshift. Thus we adopt Vexp — 200 km s~^ as a typical value. 

9. The bulk expansion is assumed, although there can be radial velocity gradient in the 
supershell. The real supershell must be complex, but we restrict the current scope of 
study within cases of bulk expansion. Here the bulk expansion implies that dust and 
gas are also strongly coupled. 
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10. We assume that dust is distributed uniformly in the ionized bubble, as well as within 
the supershell. When we let Rd be the inner radius of the dust shell, Rd can be less 
than Rmin- The outer radius of dust shell is assumed to be equal to Rmax- Although 
the dust size distribution and the number density are surely altered by the shock waves, 
we assume in this paper that they are same both in the ionized bubble and within the 
supershell. 

11. We defined a radial dust opacity through the dust shell Td-, in order to reconcile with the 
usual way of estimating dust opacity by analyzing the extinction of the UV continuum 
from central stars in the central region of Hll region. Another important physical 
parameter related with dust abundance is dust-to-gas ratio. In particular, we define 
the dust-to-Hl ratio in the neutral supershell by D = Td/ {NY{i(Td{^ — ^)}{(-Rmax ~ 
Rmin) /{Rmax ~ Rd)}^d/ ^Hh Here the typical mass of grains is defined by = 
4:7r/3a^Pd, Od is the absorption cross-section of dust grains, A is albedo, mni is the mass 
of hydrogen atom, and is the radius of dust grains effective at 1216A. We adopt the 
mass density of dust grains pd = Ig cm~^ and the size of dust grain ad = 1.9 x 10~^ cm. 
We adopt the typical albedo ^ = 0.5 in accordance with Draine & Lee (1984). 

We use the Monte Carlo method whose detailed description was presented in our pre- 
vious papers (Ahn, Lee, & Lee 2000, 2001, 2002, 2003a). In the Monte Carlo code, we 
calculate the integrated path lengths of a Lya photon before its escape only if the photon is 
passing through the dust shell. Thus the integrated dust opacity experienced by the photon 
is given by rj = Td/ro^ - Si. Here Si is a path length in units of tq along which a Lya 
photon propagates between the (i — l)-th scattering and the i-th scattering within the dust 
shell in Rd < r < Rmax- A more detailed description of the Monte Carlo code for dust- free 
supershell was given in Ahn, Lee, & Lee (2003b). 



3. Results 

We present the emergent profiles which are calculated using the Monte Carlo method in 
Figures 1 and 2. We can see that there are two peaks for each case, which is similar to our 
previous study (Ahn, Lee, & Lee 2003a,b). We call the peaks at F ~ Vexp the primary peaks, 
and the peaks a,t V W^xp the secondary peaks. Figure 1 show the emergent profiles for 

A^Hi = 10^® cm~^ and A^hi = 10^° cm~^. When we compare the left panels with Rd = OARmax 
to the right panels with Rd = O.dRmax, we can sec that the secondary peaks are destroyed 
efficiently if the dust shell extends into the ionized bubble (i.e. cases with Rd = OARmax)- 

This fact means that the emergent Lya photons, forming the peaks, have experienced 
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almost the same amount of dust extinction, if dust is distributed only within the supershell. 
This can be explained as follows. Every back-scattering causes rcdshift of Lya photons. 
Since the supershell is optically very thick, Lja photons are back-scattered at the shallow 
inner part of the supershell. Hence, the dust opacity contributed by path length within the 
supershell is relatively small, and the contribution during the final escape dominates the 
total dust opacity. Hence, the total dust opacity for most of Lya photons escaped after 
back-scattering becomes similar, which results in the inefficient decrease of secondary peaks. 

The situation becomes different when the dust shell extends into the ionized bubble, 
Rd < Rniin- While they arc scattered back and forth in the bubble by back-scattering, the 
Lya photons are redshifted and eventually escape to form the secondary and the tertiary 
peaks. Thus, back-scattered Lya photons have large path length within the dust shell in the 
ionized bubble. The dust opacity contributed by the dust shell within the ionized bubble is 
roughly proportional to the number of back-scattering. Hence the secondary and the tertiary 
peaks are more rapidly decreased than the primary peak. 

The survival fraction of Lya photons is mostly affected by the dust opacity within the 
supershell, which is determined by two factors; one is the dust opacity within the supershell, 
and the other is the integrated path length within the supershell. According to theoretical 
researches (Adams 1972; Harrington 1973; Neufcld 1990; Ahn, Lee, & Lee 2002), we usually 
call the medium extremely thick when aro > 10'^. In this case, a Lya photon experiences a 
large number of resonance core scattering, during which it happens to scatter with a hydrogen 
atom fastly moving at the tails of the Boltzmann velocity distribution. Then, the photon 
becomes a wing photon, for which the medium is transparent. The wing photon restores 
gradually to the core photon while it experiences a series of wing scattering, during which the 
photon wanders around the medium. This wandering can be regarded as a spatial diffusion, 
and can be described by random walk process. The photon has now core frequency again, 
and it repeats such processes. We call this one cycle an excursion. During such an excursion, 
Ly« can escape the medium. As a corollary, Lya photons suffer a larger amount of dust 
extinction either when the integrated path length of those photons within the supershell is 
increased or when the dust opacity within the supershell is very high. 

Firstly, the dust opacity within the supershell is rd{Rmax — Rmin)/{Rmax ~ Rd}- Thus, 
although the HI column densities are same, the dust opacity within the supershell increases as 
Rd approaches to Rmin- As a result, in Figure 3, the survival fraction curves for R^ — 0.9Rmax 
are steeper than those for Rd = OAR^ax for cases with the same A^hi and b. 

Secondly, as otq increases, both the number of wing scattering per excursions and the 
number of excursions within the supershell increases. Hence, the path length within the 
supershell is an increasing function of otq. In Figure 3, for cases with the same Doppler 
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parameter 6 = 80 km s ^, we can see that the curve with A^hi = 10^° cm ^ is decreased 
much faster than that with A^hi = 10^^ cm~^. 

Thirdly, we know in Section 2 that that a oc b"^ and Tq oc b~^. Thus, qtq oc b~^. Since 
the number of excursion within the supershell is proportional to otq. We can expect that the 
destruction of the secondary and higher peaks is enhanced for small b. We show in Figure 2 

the changes of profiles for different Doppler widths, when both the column density of neutral 
hydrogen iVni and the dust configuration (or and Td) are same. The secondary peak 
is decreased more rapidly as b is decreased. In addition, the width of the primary peaks 
decreases for small b. The red part of the primary peak is formed of the photons which are 
scattered backward. Those photons experiences dust extinction both within the supershell 
by resonance scattering and in the dust shell. Thus they can be easily destroyed by dust 
during those scattering processes, and so the primary peak becomes narrow. We show in 
Figure 3 the variation of survival fraction of Lya photons for the cases of different Doppler 
parameters, while A^hi, Td, and Rd are same. When we compare the solid line with solid 
dots with the dotted line with solid dots, we can see that the curve for cases with small b 
decreases rapidly. 

On the other hand, we can see in Figures 1 and 2 that both the widths of peaks and 

the velocity difference between peaks are largely insensitive to dust opacity. We showed that 
the kinematics of outfiowing material is imprinted on both the velocity width of the peaks 
and the velocity differences between peaks (Ahn, Lee, & Lee 2003a,b). Therefore, those 
characteristics can be used potentially to estimate the expansion velocity of the supershell 
even in dusty interstellar media. 

4. Summary 

In this paper we emphasize that the spatial distribution of dust grains in the supershell 
is as important in forming the P-Cygni type Lja emission profiles as the kinematics of the 
supershell. The escape probability of Lja photons is greatly enhanced when the supershell is 
expanding with respect to the central Lya sources. That is, the Lya resonance scattering is 
off-centered and the escape probability of Lya photons is enhanced. As the number of back- 
scattering increases, Lya photons are redshifted more and more, which enhances the escape 
probability. Lya photons are scattered backward by high HI column density of expanding 
supershell. Therefore, in the cases of dust-free supershell, these back-scattering processes 
result in a series of peaks in the redward of Lja. However, majority of observed Lya emission 
show no such peaks. In this paper, we have found that the single peak can be formed, only 
when dust extends in the ionized bubble, as well as within the supershell. In this case. 
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the integrated path length in dust shell increases as Lya photons scatter back and forth 
by the shallow part of the inner wall of the supershell. Thus, the integrated dust opacities 
increases as the inner radius of the dust shell decreases, which is main discovery of this study. 
Moreover, when Lya photons suffer a large number of resonance core scattering and wing 
scattering within the supershell, their path lengths become very large. Thus, the dust opacity 
within the supershell is also one of decisive factor that determines the overall escapability 
of Lya photons. Furthermore, both the optical depth and the Doppler parameter of the 
supershell are another crucial factor to determine the survival fraction of Lya photons. We 
have found that the kinematic information of outflowing media surrounding Lya sources is 
imprinted on the Lya proflles, and it is conserved even when the medium is dusty. However, 
our model is very simplified one, whose limitation is described in Section 2 in detail. More 
work is needed in the future. 
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Fig. 1. — Profiles of emergent Lya emission. The upper panel is for cases with iVni = 
10^° cm~^, and the lower panel is for cases with A^hi = 10^^ cm^^. The inner radii of cases 
in the left panels are OARmax, and those of left panels arc 0.9R„iax- The expansion velocity 
of the supershell Vexp = 200 km s~^, and the Doppler parameter 6 = 80 km s^^ in common. 
In order to show the relative flux ratio between the primary and the secondary peaks, all 
the profiles are normalized with respect to the primary peak. The vertical coordinate means 
fiux, and the survival fraction of Lya photons is shown in Figure 3. The dust-to-gas ratios 
in the supershell for Td — 1 are given in each box. Note that the dust-to-(HI-l-He) ratios 
of Our Galaxy and the Large Magellanic Cloud are 0.007 and 0.002, respectively (Whittet 
1992). 
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Fig. 2. — Profiles of emergent Lya emission for cases with different Doppler widths and same 
other parameters (A^m — 10^^ cm~^, Vg^p = 200 km s~^, — OAR^ax)- All the profiles are 
normalized with respect to the primary peak. The relative strength is represented by dotted 
lines of survival fraction curves in Figure 3. Note that the secondary peaks decreases more 
sensitively as b decreases. 



-13- 




Fig. 3. — Survival fraction of Lya photons with dust opacities. The model parameters of 
the supershell is given in the figure. We can see that more Lja photons are destroyed as 
Rd increases to reach Rmin- However, we see in Figures 1 and 2 that the secondary peak is 
more rapidly decreased as Rd decreases. 



